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Abstract. This paper intends to review some of the prominent properties of strongly coupled classical 
plasmas having in mind the possible link with the quark-gluon plasma created in heavy-ion collisions. 
Thermodynamic and transport properties of classical liquid-state one-component plasmas are described 
and features of collective excitations are presented. 



PACS. 52.27.Gr Strongly-coupled plasmas 
52.25.Fi Transport properties 



52.27.Lw Dusty or complex plasmas; plasma crystals 



' 1 Introduction 

^ In the RHIC experiments Au+Au collisions at ultra- 
• relativistic energies take place and an extremely high en- 
' ergy density system is created. The experiments demon- 
, strate that a collective state has been created in these 

■ nucleus-nucleus collisions, where matter consists of a large 
, number of deconfined quasi-free constituents of the nucle- 

■ ons, namely quarks and gluons. The striking discovery was 
' that these particles are, however, in a strongly interacting 
. phase, resembling a liquid rather than a gas: this phase 
' is being referred to as a strongly coupled Quark Gluon 
; Plasma (sQGP) [1]. 

The strongly coupled quark-gluon plasma is in many 
' ways similar to certain kinds of conventional (electromag- 
i netic) plasmas consisting of electrically charged particles 
' (electrons, ions or large charged mesoscopic grains), which 
also exhibit liquid or even solid-like behavior. These plas- 
mas are known as strongly coupled plasmas and are char- 
acterized by an inter-particle potential energy which dom- 
inates over the (thermal) kinetic energy of the particles. 
Strongly coupled plasmas occur in electrical discharges, 
in cryogenic traps and storage rings, in semiconductors, 
and in astrophysical systems (interior of giant planets 
and white dwarfs). Investigations of these physical sys- 
tems have been a major field of activity for some time 

Plasmas are extremely versatile, as illustrated in 
Fig. [11 which shows some plasma types over the density 
- temperature plane. Besides the more conventional types 
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Fig. 1. Different types of plasmas over the density - tempera- 
ture plane. Note the extremely wide range of these parameters. 
Strongly coupled plasmas are located right from the F = 1 line. 
{r characterizes the ratio of the potential to kinetic energy, see 
Eq. Q.) 



of plasmas, the approximate location of the sQGP is also 
shown in Fig. [T] 

Strongly-coupled plasmas may consist of different 
charged species. Neutron star crusts are composed of fully 
stripped iron ions and electrons. Besides electrons, in the 
core of Jovian planets we find a binary mixture of and 
He^+ ions [3j, while the core of white dwarf stars con- 
sists of a mixture of fully stripped ions of C, N, and O [1]. 
Dusty plasmas, in addition to electrons and ions, also con- 
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tain another component of mesoscopic dust grains, which 
charge up and respond to electromagnetic fields similarly 
as electrons and ions ,5,6,7,8,19]. 

Part of the different plasmas listed above can be 
described within the framework of the one- component 
plasma model (see e.g. [10],) which considers explicitly 
only a single type of charged species and uses a poten- 
tial that accounts for the presence and effects of other 
types of species. This latter may considered as a charge- 
neutralizing background, which is either non-polarizable 
or polarizable. In the former case the interaction of the 
main plasma constituents can be expressed by the 



(r) 
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Coulomb potential energy, whereas in the case of polariz- 
able background the use of the 







exp(-r/AD) 
47r£:n r 



(2) 



Yukawa potential energy is appropriate to account for 
screening effects {Q is the charge of the particles and Ad is 
the Debye length). As examples of systems for which the 
Yukawa potential can be used, dusty plasmas t5j6,.7j.8;S] 
and charged colloids p!Tl[T2l[T3lfT4] may be mentioned. 

Strongly coupled plasmas appear in nature and in lab- 
oratory environments in both three-dimensional (3D) and 
two-dimensional (2D) settings. While 3D systems are more 
widespread, notable examples of 2D systems are the layer 
of dust particles levitated in gaseous discharges [TMTC] 
and the layer of electrons over liquid helium surface [17[ 

[ig. 

Strongly coupled one-component Coulomb systems are 
fully characterized by the coupling parameter. 



r = 



1 



(3) 



where a is the Wigner-Seitz (WS) radius, and T is the 
temperature. In the case of Yukawa interaction an addi- 
tional essential parameter is the screening parameter. 



Ar 



(4) 



The coupling parameter F is the measure of the ratio 
of the average potential energy to the average kinetic en- 
ergy per particle. The strong coupling regime corresponds 
to r > 1. In the K — > limit the interaction reduces to 
Coulomb type, while at k ^ oo it approximates the prop- 
erties of a hard sphere interaction. 

Recent work indicates that the coupling parameter for 
the sQGP is expected to be in the order of one [19] (see 
Fig.[T]). This is exactly the reason why methods tradition- 
ally used in the mathematical description of strongly cou- 
pled plasmas may become useful in the physics of sQGP. 
Therefore, reviewing the prominent properties of more 
conventional types of strongly coupled plasmas is of in- 
terest and this is indeed the motivation of the present 



paper. We, on the other hand, do not go beyond the 
one-component plasma (OCP) model, while a possible de- 
scription of the sQGP would clearly require a multicom- 
ponent plasma model. The methods applicable to single- 
component systems (such as the ones we deal with in the 
present paper) serve as the basis of the description of mul- 
ticomponent systems, such as ionic mixtures [3]l4l[20] and 
charged particle bilayers [H] . 

At present, large enough scale numerical simulations, 
to reproduce e.g. collective excitations, are only feasible 
for classical systems. Therefore we also restrict our stud- 
ies (presented here) to classical systems. More sophisti- 
cated, but still classical models aiming partial descrip- 
tion of sQGP phenomena including electric and magnetic 
charges were developed by the group of Shuryak |22j . 
An attempt to include non-Abelian color-color interaction 
into the classical simulation was presented in [23] . It is ex- 
pected that with the advance of computational resources 
large scale simulations for quantum systems will become 
realistic within the coming decade. 

In Section 2 we introduce the theoretical and numer- 
ical methods applied in our studies. Section 3 describes 
basic thermodynamic properties of the strongly coupled 
one-component plasma (sOCP). Section 4 deals with the 
transport properties of sOGP, while Section 5 presents 
properties of collective excitations characteristic for the 
liquid phase sOGP. Finally, Section 6 gives a short sum- 
mary of the paper. 



2 Theoretical and numerical methods 

Many body systems can be treated theoretically in a 
straightforward way in the extreme limits of both weak in- 
teraction and very strong interaction. In the first case, one 
is faced with a gaseous system, or a Vlasov plasma, where 
correlation effects can be treated perturbatively (I^ <C 1). 
Sophisticated theoretical approaches, like diagrammatic 
expansions [24] make it possible to extend standard meth- 
ods to obtain thermodynamic results in the moderately 
coupled regime. The random phase approximation (RPA) 
[2S] method, based on the linear response theory, is a use- 
ful tool to calculate dynamical properties (wave disper- 
sions) in the case when correlation effects are negligible. 

In the case of very strong interaction, the systems crys- 
tallize, the particles are completely localized and phonons 
are the principal excitations. For such conditions lattice- 
summation techniques serve as solid basis to obtain wave 
dispersion information. 

In the intermediate regime - in the strongly coupled 
liquid phase - the localization of the particles in the lo- 
cal minima of the potential surface still prevails, how- 
ever due to the diffusion of the particles themselves the 
time of localization is finite [26]. A successful theoreti- 
cal approach for calculating structural properties, like the 
static structure function S{k), is the Hyper-Netted-Chain 
(HNG) method with the Percus-Yevik (PY) equation [27] . 

The localization of the particles (which may typically 
cover a period of several plasma oscillation cycles) serves 
as the basis of the Quasi- Localized Charge Approximation 
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(QLCA) method [281129], The QLCA uses structural in- 
formation (in the form of static structure function S{k) or 
the pair correlation function g{r)) as input for the calcu- 
lations of the dispersion relations of the collective modes. 
The conceptual basis for the QLCA has been a model that 
implies the following assumptions about the behavior of 
strongly coupled Coulomb or Yukawa liquids [28^: (i) in 
the potential landscape deep potential minima form that 
are capable of trapping (caging) charged particles; (ii) a 
caged charge oscillates with a frequency that is determined 
both by the local potential well and the interaction with 
the other (caged) particles in their instantaneously frozen 
positions; (iii) the potential landscape changes slowly to 
allow the charges to execute a fair number of oscillations; 
(iv) the escape from the cages of the particles is caused 
by the gradual disintegration of the caging environment; 
the timescale of this process is determined by the cou- 
pling strength F; (v) the (time and velocity dependent) 
correlation between a selected pair of particles is well ap- 
proximated by the (time and velocity independent) equi- 
librium pair correlation; (vi) the frequency spectrum cal- 
culated from the averaged (correlated) distribution of par- 
ticles represents, in a good approximation the average of 
the distribution of frequencies originating from the actual 
ensemble. These assumptions have been confirmed to be 
reasonable in a series of studies (e.g. [261129]). 

The main concern of the QLCA theory is the analy- 
sis of the collective behavior in strongly coupled many- 
particle systems. The formal tools for this are the dielec- 
tric function e^^(k, w) [having a tensor character (sub- 
scripts) in real space and a matrix character (superscripts) 
in species space] and the dynamical structure function 
S^^ {]i,uj), or more generally, the dynamical current- 
current correlation function T^^(k, w). The principal ap- 
proximation of the QLCA method is to replace the fluc- 
tuating microscopic densities and their products by their 
ensemble averages, making use of the S'(k:) static structure 
function of the system. 

All the information pertaining to the mode structure is 
contained in the dielectric matrix that has a longitudinal 
and a transverse element: 



eL/T(k,t^) = 1 - 



^o'(k) 



-i?L/T(k) 



(5) 



Here the -DL(k) and -DT(k) local field functions are the 
respective projections of the QLCA dynamical matrix 
D^^(k) ^28j . which is a functional of the equilibrium pair 
correlation function (PCF) h{r) = g{r) — 1 or its Fourier 
transform h(k): 

D,,{k) = f d2rAf^,(r)[e'k-r _ ^^ir) (6) 
m J 

with Mf^^{r) = d^d,y(f){r) being the dipole-dipole interac- 
tion potential associated with (/)(r). 

The longitudinal and transverse modes are now deter- 
mined from the dispersion relations 



eL(k,w) = 0, £T^(k,^^) 0. 



(7) 



Besides the theoretical approaches computer simula- 
tions have proven to be invaluable tools for investigations 
of strongly coupled liquids of charged particles. Monte 
Carlo (MC) and molecular dynamics (MD) methods have 
widely been applied in studies of the equilibrium and 
transport properties, as well as of dynamical effects and 
collective excitations. The main difference between the 
two techniques is that in a MC simulation independent 
particle configurations of a canonical ensemble are gener- 
ated, whereas MD simulations provide information about 
the time-dependent phase space coordinates of the parti- 
cles, this way allowing studies of dynamical properties. 

Molecular dynamics simulations follow the motion of 
particles by integrating their equations of motion while 
accounting for the pairwise interaction of the particles, see 
e.g. [3^. Assuming that the dynamics is Newtonian, for 
each of the particles the Newtonian equation of motion, 

mv, = F, (8) 

has to be solved. Here is the total force acting on the 
i-th particle due to all the other particles and due to any 
external (e.g. electric and/or magnetic) field. The way how 
Fi is calculated will be discussed below. In some physical 
systems further details have to be considered beyond the 
Newtonian approximation. As an example charged colloids 
can be mentioned where the Brownian molecular dynam- 
ics simulation [121l31j is widely used. Brownian molecular 
dynamics simulations take into account in the equation of 
motion solvent friction and a random Langevin force R(t) 
acting on the particles. 

In the rest of the paper we restrict our studies to sys- 
tems where Newtonian dynamics is a reasonable approx- 
imation. In the case of short-range potentials the calcu- 
lation of the force acting on a particle of the system, F^, 
is relatively simple. In this case MD methods make use of 
the truncation of the interaction potential thereby limiting 
the need for the summation of pairwise interaction around 
a test particle to a region of finite size. In the case of long- 
range interactions (e.g. Coulomb or 1ow-k Yukawa poten- 
tials), which are of interest here, however, such truncation 
of the potential is not allowed, and thus special techniques, 
like Ewald summation [32], have to be used in MD simu- 
lations. Besides the Ewald summation technique there ex- 
ist few additional methods, like the fast multipole method 
and the particle-particle particle-mesh method (PPPM, or 
P3M), which can be used to handle long-range interaction 
potentials, see e.g. [33]. It is this latter - widely used [331 
l35] - method, which we also choose as the simulation ap- 
proach for our studies presented here. The PPPM method 
was originally applied for the case Coulombic interaction 
"36]. In the PPPM scheme the interparticle force is parti- 
tioned into (i) a force component -FpM that can be calcu- 
lated on a mesh (the "mesh force" ) and (ii) a short-range 
("correction") force Jpp, which is to be applied to closely 
separated pairs of particles only. In the mesh part of the 
calculation charged clouds are used instead of point-like 
particles and their interaction is calculated on a compu- 
tational mesh, taking also into account periodic images 
(for more details see [36]). This way the PPPM method 
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makes it possible to take into account periodic images of 
the system (in the PM part), without truncating the long 
range Coulomb or 1ow-k Yukawa potentials. (For high k 
values the PP part alone provides sufhcient accuracy, in 
these cases the mesh part of the calculation is not used.) 



3 One-component plasma properties 

Starting with the pioneering work of Brush, Sahlin and 
Teller jlO] and followed by the systematic studies of 
Hansen and coworkers [37 ( l38 l [39 l l40] properties of one- 
component plasmas have been explored by computer sim- 
ulation and theoretical approaches. 

In three dimensions (3D) at k = the liquid phase 
is limited to coupling parameter values T < 175 [H]. A 
first order phase transition was identified to take place at 
7^ ~ 175, where the plasma was found to crystallize into 
hcc lattice [H]. We note that at k > 0, the 3D systems 
may crystallize either in bcc or in fee lattices, depending on 
the value of k |43j . Crystallization of the plasma was also 
experimentally confirmed to take place in many different 
systems, e.g. space plasmas [l^, in laser-cooled trapped 
ion plasmas [l^ and expanding neutral plasmas [46j, as 
well as in ion storage rings [47ll48] . 

In two dimensions (2D) crystallization into hexagonal 
lattice occurs at a lower value of coupling, at « 137, as 
found by both computer simulations [49] and by experi- 
ments jl7l . In two dimensions the crystallized form of the 
systems is always hexagonal. While most of the studies on 
2D system have been carried out in the crystalline state 
( "plasma crystals" ) the liquid state also receives more at- 
tention nowadays [5CTII5T] . 




Fig. 2. Pair correlation functions of the classical 3D OCP. 
Dependence on F in the Coulomb case (a), and dependence 
on the screening parameter of the Yukawa potential at fixed 
coupling r — 160 (b). 



3.1 Liquid state properties 

At high values of the coupling coefficient plasmas exhibit 
strong structural correlations. Such correlations can easily 
been studied by examining the pair correlation function 

(9) 

of the system {V is the volume of the system consisting of 
TV particles). Figure [2lja) shows pair correlation functions 
for the Coulomb OCP for a series of coupling parameter 
values, while Fig. Wih) illustrates the changes of g{r) as a 
consequence of screening (k > 0). With decreasing F the 
peak amplitudes of the pair correlation function decrease, 
but the positions of the peaks remain nearly unchanged. 
This remarkable feature of the pair correlation functions 
indicates that the local environment of the particles in 
the liquid phase still resembles the underlying {F — > oo) 
lattice configuration. Increasing screening also decreases 
the degree of correlation as it can be seen in Fig. [Hb). 

Particle simulation methods provide the full informa- 
tion about the mechanical state of the system (position 
and velocity of every particle in every discrete time step) . 



However, to extract macroscopic quantities (like thermo- 
dynamic energy, pressure, compressibility) it is useful to 
construct first distribution functions (like the g[r) pair 
correlation function, being specially important for liquid- 
state studies, where overall isotropy is a natural assump- 
tion) and use the standard tools of classical statistical me- 
chanics (e.g. [S2]) to calculate thermodynamic character- 
istics. 

The total energy of the system can be expressed as: 
E^^Nk^T + U, (10) 

which results in 

^ = 2^Br+ I V{r)9{r) Anr^dr, (11) 

where U is the excess energy, n is the particle number den- 
sity and '^(r) is the interaction pair-potential. The pres- 
sure can be calculated using 

p = nk^T r ^.g(r) Anr^dr, (12) 

6 Jo or 
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and the isothermal compressibihty is expressed by 



kv.T 



dn 
dp 



[g{r) - 1] Airr'^dr. (13) 



The apphcation of the above formulae, derived from 
fundamental principles of statistical mechanics, depends 
on the form of the (p{r) pair potential. For Coulomb inter- 
action {(p{r) ~ 1/r) the integral in Eq. (fTT|) is divergent, if 
the contribution of the uniform neutralizing background 
is not properly taken into account. This can be done sim- 
ply by replacing g{r) with h(r) = g(r) — 1. For Yukawa 
interaction the contribution of the polarized background 
is finite {Eh ~ 1/k^, 3D Hartree energy). The potential 
energy can be written as the sum of the Hartree energy 
and the correlational energy: 



^Yukawa 

N 



4:TT£q a 



\+ I ft.(f)fe~'"'df 



(14) 



where r — r/a. 



4 Transport properties 
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Fig. 3. Self-diffusion coefficient of the 3D classical one- 
component plasma (OCP). Hansen et at |39j . OH: Ohta and 
Hamaguchi [53], DKG: Donko, Kalman and Golden, calcu- 
lated from cage correlation functions [26], DHK: present MD 
data. The self-diffusion coefficient has been normalized as: 
D* = D/a^Loo- All data correspond to k = 0, except OH, 
which is for n = 0.1. 



In this section we review the data available in the lit- 
erature for the basic transport coefficients (self-diffusion, 
shear viscosity and thermal conductivity) of the classical 
one-component plasma. We mainly present data here for 
Coulomb OCP, although some data obtained for Yukawa 
potential at very low values of the screening parameter 
(0.01 < K < 0.1) are also shown for comparison. Trans- 
port properties of Yukawa systems characterized by such 
low K values are very close to those of Coulomb systems. 
As already mentioned in the introduction, there is a grow- 
ing interest in Yukawa systems. (For more details about 
transport properties of Yukawa systems the Reader is re- 
ferred to the original publications.) 

Simulation techniques have become indispensable tools 
for the determination of transport coefficients. The two 
main approaches of molecular simulation are the equi- 
librium and the non- equilibrium methods. In the former 
one the transport coefficients are derived from correlation 
functions of microscopic quantities using the Green-Kubo 
(GK) relations. In non-equilibrium simulations a pertur- 
bation is apphed to the system and the system's response 
is measured. 



4.1 Transport coefficients of the classical 
one-component plasma 

4.1.1 Self-diffusion 

In equilibrium molecular dynamics simulations one can 
compute the self-diffusion coefficient D from cither the 
Green-Kubo relation 

^ = ^^"(v(Ov(0))dt, (15) 



(i.e. via the velocity autocorrelation function), or from the 
Einstein formula: 



D= hm -(|r,(i)-r,(0)p}, (16) 
t— >cx) or 

(i.e. from the mean square displacement of the particles). 
In the above formulae averaging is taken over particles 
and different initial times. 

The known for the self-diffusion coefficient are dis- 
played in Fig. [21 The data have been normalized according 
to D* — D/a^ujQ, where loq = \/ riQ"^ / e^m is the plasma 
frequency and m is the mass of theparticles. 

Hansen et al. made use of (|15l) to obtain the self- 
diffusion coefficient of the Coulomb OCP. Their results 
were found to follow the approximate relation D* = 
2.95r-i-3'' [M]- Ohta and Hamaguchi [53] obtained the 
self- diffusion coefficient for Yukawa liquids from molecu- 
lar dynamics simulations using |[T6)) . Their results for k = 
0.1 as well as our present data (based on the same com- 
putational procedure) obtained for k — are also shown 
in Fig. [3] These more recent MD data fall very close to 
the those given by the above formula. 

An additional set of data derived on the basis of the 
caged behavior and jumping of the particles in the strongly 
coupled liquid phase [26] is also shown in Fig. [3] This data 
set agrees surprisingly well with the results of the "direct" 
MD calculations for D over a wide domain of the coupling 
parameter. 
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4.1.2 Shear viscosity 

The shear viscosity data for the 3D OCP are shown in 
Fig. m The shear viscosity coefficient has been normahzed 
as: r]* = rj/mna'^ujo. 

The shear {ij) and bulk {() viscosity of the 3D OCP 
was first derived by Vieillefosse and Hansen )40J from the 
transverse and longitudinal current correlation functions 
of the plasma. They have found that the shear viscosity 
exhibits a minimum at w 20. The other main finding of 
their work was that the bulk viscosity is orders of mag- 
nitude smaller compared to the shear viscosity. The cal- 
culations of Wallenborn and Baus [5HI55] were based on 
the kinetic theory of the OCP to calculate 77. Their results 
were in a factor of three agreement with the previous re- 
sults [30] at r = 1 and within a factor of two agreement at 
r = 160. The minimum value of rj agreed well for both re- 
ports, however the position of the minimum was reported 
in [54] to occur at a lower coupling value, F 8. Molecular 
dynamics simulation was first applied by Bernu, Vieille- 
fosse and Hansen [5S1I57] to obtain transport parameters 
through the Green-Kubo relations. Donko and Nyi'ri [58] 
used a non-equilibrium MD simulation technique to deter- 
mine the shear viscosity, while subsequently, Bastea [59j 
applied equilibrium simulation and obtained rj from the 
Green-Kubo relation. Daligault [60j has found that the 
shear viscosity of the OCP follows an Arrhenius type be- 
havior at high r values. This is shown in Fig. [4] by dashed 
line. 

Salin and Caillol |61| have carried out equilibrium 
molecular dynamics computations for the shear and bulk 
viscosity coefficients, as well as for the thermal conduc- 
tivity of the Yukawa one-component plasmas. They have 
implemented Ewald sums for the potentials, the forces, 
and for all the currents which enter the Kubo formulas. 
Saigo and Hamaguchi 62J have also used the Green-Kubo 
relations for the calculations of the shear viscosity. As a 
refinement, the effect of the plasma environment is dusty 
plasmas has been taken into account through Langevin 
dynamics in the calculation of shear viscosity of Yukawa 
systems j63] . 

4.1.3 Thermal conductivity 

Thermal conductivity data for the 3D OCP are shown 
in Fig. [5l The thermal conductivity coefficient has been 
normalized as: A* = X/nkBO^LUQ. 

We have reported non-equilibrium molecular dynam- 
ics calculation of the thermal conductivity of the classi- 
cal OCP in [64] . In contrast with the studies of Bernu et 
al [56p57j where the transport coefficients were obtained 
from the simulation of an equilibrium system, we applied 
a perturbation to the system and deduced A* from the re- 
laxation time of the system towards the equilibrium state. 
Donko and Hartmann [65j applied the non-equilibrium 
MD method proposed by Miiller-Plathe [66] to calculate 
the thermal conductivity of Yukawa liquids. 

As regards to Yukawa systems, transport parameters 
have been studied in several papers. Besides the work 
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Fig. 4. Shear viscosity coefficient of the 3D classical one- 
component plasma (OCP). DN: Donko and Nyiri [58] using 
1024 and 8192 particles, WB: Wallenborn and Baus [Rll^ , 
VH: Vieillefosse and Hansen [40], BV : Bernu et al. [56][57] . 
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SH : Saigo and Hamaguchi ^2]. (The results of [60] have been 
scaled to match the minimum value of 77.) 
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Fig. 5. Thermal conductivity coefficient of the 3D classical 
one-component plasma (OCP). BVH: Bernu et al. |56II57] , DN: 
Donko and Nyi'ri [58] using 1024 and 8192 particles, SC: Salin 
and Caillol JjT, DH: Donko and Hartmann [65] . 



of Salin and Caillol [ST] (mentioned above), Faussurier 
and Murillo obtained thermal conductivity (as well as 
self-diffusion and shear viscosity) values for the Yukawa 
OCP through its mapping with the Coulomb OCP sys- 
tem, based on the Gibbs-Bogolyubov inequality ^67] . 

5 Collective behavior 

Collective excitations (waves) are prominent features of 
plasmas. Depending on the dimensionality and the con- 
finement of the system, different collective excitations 
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(longitudinal and transverse modes) may show up. Longi- 
tudinal modes can fully be characterized by the dynam- 
ical structure function S{k,Lu), while transverse modes 
can be studied through the analysis of the transverse cur- 
rent fluctuation spectra T(k,uj). The corresponding cur- 
rent fluctuation spectra for the longitudinal mode, L{k, uj), 
is linked with the dynamical structure function through 
L{k,uj) = {uj"^ /k'^)S{k,uj). Collective excitations are iden- 
tified as peaks in these spectra and dispersion relation are 
derived by observing the change of the frequency (where 
the peaks are found) with wave number. Additionally, 
the widths of the peaks in the spectra convey informa- 
tion about the lifetime of excitations (associated with the 
damping of the waves), as well as about the distribution of 
the mode frequencies due to the disordered particle con- 
figuration in the liquid phase. 

Collective effects in Coulomb [7T1I72] and Yukawa [73l 
[741175] plasmas have extensively been investigated. In 
Coulomb systems the longitudinal (plasmon) mode is 
known to have a frequency w — > Wp at fc — » 0. The trans- 
verse mode exhibits an acoustic dispersion. Rosenberg and 
Kalman [75| investigated the dispersion relation for dust 
acoustic waves in a strongly coupled dusty plasma with 
the aid of the QLCA scheme, generalized to take into ac- 
count electron and/or ion screening of the dust grains. 
Hamaguchi and Ohta [731174] studied the wave dispersion 
relations in the fluid phase of Yukawa systems through 
molecular dynamics simulations. They have demonstrated 
that the transverse wave dispersion has a cutoff at a long 
wavelength even in the case of weak screening. Their re- 
sults have confirmed the earlier theoretical predictions 
[75] . The QLCA method has subsequently been applied to 
determine the properties of the transverse (shear) mode 
in strongly coupled dusty plasmas [76] . For this mode the 
dispersion was found be characterized by a low-fc acoustic 
behavior and by a frequency maximum lying well below 
the plasma frequency. The dispersion curves of the longitu- 
dinal and transverse modes were demonstrated to merge 
at high wave number around the Einstein frequency of 
localized oscillations. Experimental observation of trans- 
verse shear waves in the strongly coupled liquid phase of 
a three-dimensional (layered) dusty plasma have been re- 
ported by Pramanik et al. Tf]. The collective modes of 
dusty plasmas in the liquid phase have also been investi- 
gated theoretically by Murillo [78 1 179 ( 180]. 



In the following we present MD simulation results for 
the collective excitations in 3D Yukawa liquids, and com- 
pare these with the predictions of the QLCA theory. The 
simulations have been carried out using N = 12800 par- 
ticles. In the MD simulation information about the (ther- 
mally excited) collective modes and their dispersion is ob- 
tained from the Fourier analysis of the correlation spectra 
of the density fluctuations 



s 
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Fig. 6. Longitudinal current fluctuation spectra in 3D 
Coulomb OCP at r = 160 and in a Yukawa OCP at T = 
200, K = 1. k = ka denotes the dimensionless wave number, 
its values are given in (b), while in (a) the arrow indicates 
increasing values of k. 



yielding the dynamical structure function as [39] : 



S{k,uj) 



1 



2ttN at^oo at 



(18) 



where AT is the length of data recording period and 
p{k,uj) = T[p{k,t)\ is the Fourier transform of ([TT]) . 

Similarly, the spectra of the longitudinal and trans- 
verse current fluctuations, L{k,uj) and T{k,u:), respec- 
tively, can be obtained from Fourier analysis of the mi- 
croscopic quantities 



N 



A(fc, t) — k Vjx (t) exp [ikxj (t)] , 

N 

T{k,t) — k Vjy (t) exp \ikxj (t)] , (19) 



N 



p{k,t) = k'y^^exp\ikxj{t)] 



(17) 



where Xj and Vj are the position and velocity of the j- 
th particle. Here we assume that k is directed along the 
X axis (the system is isotropic) and accordingly omit the 
vector notation of the wave number. The way described 
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Fig. 7. Transverse current fluctuation spectra in 3D Coulomb 
OCP at r = 160 and in a Yukawa OCP at T = 200, k = 
1. The arrows indicates increasing values of the dimensionless 
wave number k (for values see Fig. [6]). 



above for the derivation of the spectra provides informa- 
tion for a series of wave numbers, which are multiples of 
fcmin = 27r/if, where H is the edge length of the simula- 
tion box. The collective modes are identified as peaks in 
the fluctuation spectra. The widths of the peaks provide 
additional information about the lifetimes of the excita- 
tions: narrow peaks correspond to longer lifetimes, while 
broad features are signals for short lived excitations. 

Representative longitudinal and transverse current 
fluctuation spectra, L{k,uj) and T{k,Lu), respectively, are 
plotted in Figs. [S] and [7| for wave numbers, which are mul- 
tiples of fcmin = fcrninO. = 0.167. L{k,ijj) obtained for the 
Coulomb case {F — 160, n — Q) peaks very nearly at the 
plasma frequency ujq . In the presence of screening (Yukawa 
potential), as shown in Fig. [6l^b), the behavior of L{k,uj) 
changes significantly: at fc — » the wave frequency lo —>■ 
0. The contrast between the k = and the k > cases 
is also well seen in Fig. |8l where the dispersion curves de- 
rived from the fluctuation spectra are displayed. The dis- 
persion curves for k > are quasi-acoustic {uj/ujq oc A:^/^), 
with a linear portion near fc = 0, which gradually extends 
when K, is increased. The {r,K) pairs for which the dis- 
persion graphs are plotted in Fig. [5] have been selected to 
represent a constant "effective" coupling 7^* = 160. This 




Fig. 8. Dispersion relations for the (a) longitudinal and (b) 
transverse modes of 3D Coulomb and Yukawa plasma liquids. 
Symbols represent molecular dynamics results, while the lines 
correspond to the predictions of the QLCA theory. The {F, ft) 
pairs are given in the legend of panel (b) . 



deflnition of F* relies on the constancy of the flrst peak 
amplitude of the pair correlation function g{f), similarly 
to the case of 2D Yukawa liquids [81] . 

Compared to those characterizing the C mode, peaks 
in the T mode spectra are rather broad, as it can be seen 
in Fig. [Tja) and (b), for the Coulomb and Yukawa cases, 
respectively. In the case of this mode there is no significant 
change between the behavior when k changes from zero to 
a nonzero value, only the mode frequency decreases, as can 
be observed in Fig.[5]^b). 

Comparison of the dispersion relations obtained from 
the MD and QLCA results [76] is presented in Fig. [8] 
The QLCA equations for the mode frequencies need the 
pair correlation function as input data. The data shown in 
Fig. [5] were obtained using MD-generated g{r) functions. 
The agreement between the (MD and QLCA) dispersion 
curves is excellent for the £ mode, while some differences 
in the frequency of the T waves can be seen in Fig. |5l^b). 
This latter may originate from the inaccurate determina- 
tion of the peak positions of the rather spread T{k,uj) 
spectra. Another difference is the cutoff of the T mode 
dispersion curve at finite wave numbers. This disappear- 
ance of the shear modes for A: — > is a well known feature 
of the liquid state [5M5^1l55j while the sharp cut-off oj — > 
for a finite k has also been observed in simulations of 
Yukawa systems [73][80] . It is noted that this cutoff is not 
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accounted for by the QLCA, as it does not include damp- 
ing effects. 

6 Summary 

This paper intended to review some of the important prop- 
erties of strongly coupled plasmas (within the framework 
of the one-component plasma (OCP) model), which might 
be relevant to the studies of strongly interacting quark- 
gluon plasma (sQGP). 

At high values of the coupling coefficient (-T 1) the 
one-component plasma exhibits liquid state properties. In 
this domain the pronounced peaks of the pair correlation 
function indicate the presence of strong correlation effects. 
The transport coefficients - self-diffusion, shear viscosity 
and thermal conductivity - have thoroughly been investi- 
gated and their behavior is quite well understood. While 
the self-diffusion coefficient decreases monotonically with 
increasing the shear viscosity and thermal conductiv- 
ity exhibit a minimum at moderate values of coupling 
(10 < r < 20), due to the different temperature depen- 
dence of the kinetic and potential contributions to these 
transport coefficients. Two types of collective excitations 
- a longitudinal and a transverse mode - have been iden- 
tified in the OCP system, the emergence of latter of them 
being attributed to strong correlations. 
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